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The factorization problem of ^-exponential distribution within nonexten- 

sive statistical mechanics is discussed on the basis of Abe's general pseudoad- 
ditivity for equilibrium systems. It is argued that the factorization of com- 
pound probability into product of the probabilities of subsystems is nothing 
but the consequence of existence of thermodynamic equilibrium in the inter- 
acting systems having Tsallis entropy. So the factorization does not needs 
independent noninteracting systems and should be respected in all exact cal- 
culations concerning interacting nonextensive subsystems. This consideration 
makes it legitimate to use ^-exponential distribution either for composite sys- 
tem or for single body in many-body systems. Some known results of ideal 

gases obtained with additive energy are reviewed. 
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I. INTRODUCTION 

In thermostatistics, the factorization of the compound probability into product of single 
body probabilities 

N 

P = II Pn, (!) 

n=l 

where p is density operator and N is the number of bodies in the system, is crucial for 
applications of statistical mechanics to many-body systems and for the statistical interpre- 
tation of thermodynamics. Boltzmann-Gibbs statistics (BGS), in considering only short 



range interactions, makes this factorization a natural result of its exponential distribution 
with additive energy. In nonextensive statistical mechanics (NSM) intended originally 
to describe complex systems with long range interactions or fractal structure of space-time, 
this factorization is not as evident as in BGS. From the beginning of this nonextensive 
theory, in order to show the nonextensive character of NSM, one has supposed Eq. (|I|) for 
a system containing statistically independent subsystems and obtained for the total entropy 
S q 1 : 

N 

ln[l + (1 - q)S q ) = J2 Ml + (1 " q)S q {n)] (2) 

n=l 

where q is the parameter of Tsallis entropy S q = — Try^ (Boltzmann constant kn — 1) 
[for iV = 2, S q = S q (l) + S q (2) + (1 — q)S q (l)S q (2) as one often finds in the literature]. As 
indicated by Tsallis [0, Eq.(|2|) expresses in fact the additivity of Renyi entropy S R = — -j^- 
if Eq.(jl]) applies. Due to this supposed "independence" of subsystems for NSM, it has 
been believed by many that exact calculations should use the additive hamiltonian H 
given by 

TV 
n=l 

where H n is the hamiltonian of n th subsystem. However, this additive hamiltonian is com- 
patible with neither Eq.(|l]) nor Eq.([|) since these relations applied to g-exponential distri- 
bution (qed) p oc [1 — (1 — q)(3H]~, as given by maximization of Tsallis entropy under some 
constraints §00, imply 

N N N k 

# = E#n + £[(<?- 1)/?]*- 1 E ( 4 ) 

n=l k=2 ni<n2<...<n fe j=l 

— H + H c , 

where (3 is the inverse temperature. So in order to keep simultaneously Eq.([|) and Eq.(|T]) 
and to apply NSM to many-body systems using one-body qed, a so called factorization 
approximation is proposed [EJ with the assumption that the second term on the right hand 
side of Eq.(f|) may be neglected. This approximation has been, explicitly or not, employed in 



the applications of NSM to the cases like, among others, quantum particle systems p^-jMj, 
turbulent flows the polytropic model of galaxies, solar neutrinos, peculiar velocity 



of galaxy clusters, electron plasma (for updated comments on some of these works, we refer 
to reference ]nj] and the references there-in). In these calculations, it is admitted that one- 
body qed applies. Although these applications clearly shows the usefulness and necessity 
of one-body qed, the approximation fl2| neglecting the correlation energy with sometimes 



weak interacting dilute particles [17| is, on the contrary, not a reassuring basis. Recently, 



it has been argued [plj that the correlation energy (H c ) given by the second term of Eq. ( 



is in general not negligible. A numerical result for iV-oscillator system || shows that the 
partition function given by using Eq.(^) is completely different from that given by using 
Eq.([|) when N is large. We are facing with a difficult question : which one of Eq.(|l]) and 
Eq. (^|) should be related to independence of subsystems? In a recent work, Beck proposed a 
nontrivial idea to define an "independence" according to Eq. (f|) [17| so we can write Eq. ([!]) 



and one-particle qed without any approximation. This means that we address correlated 
"independent systems" satisfying Eq.(J|). This idea surely needs to be justified. 

In this letter, the probability factorization problem of NSM is discussed from a different 
consideration inspired by the general pseudoadditivity of entropy |^U| and energy |^TJ re- 



quired by the existence of thermodynamic equilibrium. It will be argued that Eq.(|l]) needs 
not independence of noninteracting systems and is to be considered as a basic assumption of 
NSM. Eq.(^|) is only a kind of extensive approximation and should be employed carefully. 
On this basis, some basic applications of NSM to classic and quantum gases are revisited. 

II. THERMAL EQUILIBRIUM AND FACTORIZATION OF COMPOUND 

PROBABILITY 

Apart from one-particle qed, another issue tightly related to probability factorization is 
the establishment of zeroth law and the definition of temperature for NSM. This is obviously 
of central importance for the theory. Eq.(^) and noninteracting model due to Eq.(|I|) are 
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so deeply accepted that it was even believed that the zeroth law of thermodynamics was 
absent within NSM [[Tj| , since Eq. (|2|) does not hold with additive hamiltonian, and without 
Eq.(0), we can not talk about the zero law and temperature! Recently, a series of works have 
been published on this issue [^9j claiming the establishment of zeroth law and the definition 
of a generalized temperature on the basis of additive hamiltonian Eq.(|3|) and, I stress it, 
Eq.(^|). If this zeroth law is established by neglecting the correlation energy H c , it is merely 
approximate. If it is exact for noninteracting systems, then the paradox between Eq.(|^) and 
Eq.(P [thus Eq.(|)] persists due to qed. In addition, other questions arise : if the zeroth 
law holds only for, e.g. noninteracting particle systems, why do we have to discard BGS1 
And what is the origin of the nonextensivity? 

This confused situation and the above questions are, in our opinion, simply due to the 
fact that Eq.(|l[) is not clearly founded for NSM. It certainly implies independence of 
noninteracting systems for BGS. But does it mean the same thing for NSM1 

Very recently, Abe [^(J proposed a general pseudoadditivity for entropy required by the 
existence of thermal equilibrium in composite nonextensive systems : f(S) = f(Si) + f(S 2 ) + 
Xf(Si)f(S 2 ) where / is certain differentiable function satisfying /(0) = and A a constant 
depending on the nature of the system of interest. So for a system containing N subsystems, 
the thermal equilibrium requires following pseudoadditivity : 

N 

m[l+A/(SO] = £m[l + A/(S;)]. (5) 

71=1 

On the other hand, Eq.(|5]) applied to Tsallis entropy means f(S) = S and A = 1 — q p0| . 
which directly leads to InTrp 9 = YliLi hiTrpf or Eq.(P (i.e. (piPj) q = pjj means piPj = pij 
where pi is the probability of state i). So this probability factorization must be regarded as 
a condition or a consequence of thermodynamic equilibrium, instead of statistical indepen- 
dence. It should be raised to the rank of basic assumption for equilibrium thermodynamics 
with Tsallis entropy and must be rigorously respected by all exact calculations. In this way, 
the zeroth law becomes evident and a temperature can be straightforwardly defined with 
maximum entropy and minimum energy P, [To| , p2| . According to above discussions, all cal- 



dilations based on Eq.(P or using one-body qed matter of fact, exact applications 

of NSM. And all calculations based on additive energy Eq.(||) should now be considered as 
a kind of "extensive approximation" and should be employed with great care. 

Energy has been proved |H] to satisfy the same kind of pseudoadditivity at equilibrium 
as Eq.(|5|). If we choose f(H) = H and A = (q — l)/3, we get 

N 

ln[l + (q- l)(3H] = £ ln[l + (q - l)(3H n ] (6) 

n=l 

which is just Eq.(||) satisfying Eq.(|]). 

Here I would like to mention that, Eq.(^) or Eq.(|4]) contains inverse temperature (3. This 
surely has something to do with the fact that this relation is subject to thermal equilibrium. 
The dependence of the correlation energy H c on temperature is something of nontrivial and 
needs explanation. A possible argument is that H is only an effective Hamiltonian so the 
effects of interactions or correlations represented by q perhaps depends on the state of the 
system and thus on temperature. If this is true, we may expect T-dependent q values. 

In the following section, I will discuss some interesting consequences of Eq.(|5]). 

III. ADDITIVE Q-DEFORMED ENERGY 

Let us begin by introducing an additive "energy" for NSM for the sake of convenience. 
We refer to 

_ ln(l + (g - 1)0H) 

(S-I)/? [) 

as g-deformed hamiltonian and recast Eq.(^) into 

N 

h=Yl K- (8) 



n=l 



This means following transformations : 



e (q-l)Ph _ l e (q-l)ph n _ l 



and 
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pcx[l + (q- l^H] 1 '^ = e-P h = J] e _/3/ln . (10) 

n 

It is interesting to see that the qed can be transformed into Boltzmann-Gibbs exponential 
distribution with a deformed energy. The reader may wonder why we refer to the additive 
hamiltonian h or h n , instead of the nonadditive one H or H n , as deformed hamiltonian. It 
should be noticed that, when addressing a system of N particles, we have to write H n = 
|^ + V n for single particle so that h n = hl ^ 1+< - q ^^ n+Vn ^ ■ It is clear that H n , instead of h n , 
is the physical energy. According to this consideration, the total hamiltonian of a system of 
"ideal gas" (free particles in the sense that we do not write the correlation energy between the 
particles in Hamiltonian and let it be represented by q ^ 1 in the deformed energy) should 
be written as 

it _ _ _ _ \ Vn 1 IT /I 1 \ 

(q-l)(3 (g-l)/3 (q-l)(3 ti 2m ° 



When q = 1 (H c = 0), we recover # = E«=i £• 



I would like to emphasize here that, as shown in Eq.(|T0D, the deformed hamiltonian h 
may be employed, if we want, to recover the mathematical structure of BGS in introducing 
a g-deformed information measure I q = — In p and an entropy S = — Trpln p where p = 
(Z = Tre~P h ). We straightforwardly obtain |^ = (3 = 1/T and a deformed first law 
du = dw + TdS, where u = Trph is the deformed internal energy and dw the work on the 
system. 

IV. REVIEW OF SOME RESULTS ABOUT IDEAL GASES 

Now I will show some consequences of Eq. (0) on ideal gases. 

A. Classical ideal gas 

For classical ideal gas, the total hamiltonian should be given by Eq.(|ll|). So the total 
partition function Z is given by [in the formalism of complete distribution (Trp = 1) with 
unnormalized average energy U = TrHp q ] : 

6 



Z = Tr[l-(l-q)(3H]— (12) 

TV 

= l[Tr[l-{l-q)PH n )—< 

n=l 

= {Tx[l-(l-q)PH n ]^} N = z N 

where z = Tr[l — (1 — q)H n ]~ is the single particle partition function. The calculation of 
Z in reference p,|J can be adopted for single particle. For q > 1, we obtain : 

_ h* [q-l 2) g-d/2 _ Afi-d/2 ( -,o\ 

where d is the dimension of volume V. The total internal energy is 

d Z]- q - 1 

u= -ap^r (14) 

d - 1 



d/3 1-q 
d - 1 



dp l-q 

dN / kT \dN(l-q)/2+l _ 



2A/v( 9 -i) 
The specific heat is so given by 

_ dU _ kdN (JT)d N(i- q )/2 /.en 

We notice that the T-dependence of C v is the same as given by the extensive approximation 
but the q— and iV-dependences are changed [|[. This point is also noticed in the case of 
iV-harmonic-oscillator system ||. Another interesting point is that, in the present case, the 
variation interval of q is enlarged tol<g<l + | from 1 < q < 1 + given by the 
extensive approximation 

For q < 1, similar discussions can be made with the extension of Hilhorst formula given 
by Prato M. 



B. Quantum gas 



For nonextensive quantum gas, from the qed for grand canonical ensemble 
complete distribution and normalized average (U = Tr Hp), we can write : 

7 



,dJ with 



N 

[l + {l-q)P{H-iiN)]7=* = H[l + (l-q)/3(H n -n)]*z. (16) 

n=l 



Let this equation be replaced by 



JV 
n=l 



W = JJ e -^< (17) 



where /i' = Hi±Og»z^0] and ^ = in[i+(i-^Hn-M)] _ The total partition function £ is 



then given by 



Z = Tr[l-(g-l)/3(#-//JV)]5=* (18) 
= Tre"^' 



where e' k is the eigenvalue of h' n and the occupation number of the state k of single 
particle. For boson and fermion, we obtain, respectively, 

OO 1 1 

z=nE e ' nki3e>k = n — ^ and z = n e e_nfeK = ik 1 + e_/H )- ( i9 ) 

k n k =0 k 1 — e 6fc fc n fc =0 fc 

Then, it is straightforward to get : 

ld(\nZ) 1 1 

n x = Irpm = ---±—± = = j 20 

13 de[ e^±l [l + (l-g)/3( ei - M )]i=?±l 

where e; is the eigenvalue of the one particle hamiltonian H n . "+" and "-" correspond to 

fermions and bosons, respectively. These are just the standard quantum distributions given 



by Buyiikkilig et al with the so called factorization approximation [O and unnormalized 
average. 

On the other hand, in the present formalism of NSM, unnormalized average does not 
give simple quantum distribution similar to Eq.(pO|) for standard bosons and fermions. This 
is due to the fact that Z = U k J2n k e~ nk(3 < ^ U k [J2n k e- grikf3 <} 1/q . Here I only show that, in 
this case, Eq.([20|) becomes 

fii = lrp q n l = ^ - q (21) 

[l + (g-l)/9( C| -Ai)]5=i±l 

8 



where Q = Trp q can be regarded as a parameter depending on q. Q > 1 Q = 1 and Q < 1 
for g < 1, q — 1 and q > 1, respectively. These distributions seems interesting because they 
allow intermediate occupation number between that of bosons and fermions. In particular, 
at absolute zero, for "fermionlike" particle with "+", fix — Q when ei < u and n\ = 
when e; > /i. This means that it is possible for several "fermions" to occupy an one-particle 
quantum state if Q > 1 or q < 1. Consequently, the Fermi surface ei?atT = 0or/3 = oo 
changes as a function of Q or of the interaction between the particles : ep = ^3 where 
Ef is the conventional Fermi energy at q = 1 or Q = 1. This result can be compared to 
that of the fractional exclusion statistics (FES) |23| for intermediate particles different from 
bosons and fermions. It is not surprising to see that nonextensive statistics has similar effect 
to that of FES describing low dimension quasi-particles or elementary excitations p3| , p4| , 
because qed is just a result of interactions which are perhaps the origin of the quasi-particles. 
However, the fact that only the NSM formalism with unnormalized average can give the 
intermediate quantum distributions seems to deserve further investigation. Quite interesting 
efforts have been made by some authors ]I3, |25pE | to relate nonextensive statistics and the 
quantum distributions given by Eq.pO) to FES. 



V. CONCLUSION 

Summing up, it is argued that the nonextensive thermostatistics should be based on the 
factorization of compound probability suggested, not by "independence" of noninteracting 
systems, but by Abe's pseudoadditivity for equilibrium interacting systems having Tsallis 
entropy. So Eq.(|l]) should be viewed as a fundamental hypothesis of NSM and has to be 
rigorously satisfied by all exact calculations relative to equilibrium systems. All calculations 
based on additive hamiltonian Eq.(|3|) are to be considered as extensive energy approximation 
which should be employed with great care. In this formalism, the applications of NSM to 
ideal gases are revisited. The results are different from those given by extensive approx- 
imation. The standard quantum distributions can be given by using normalized average 
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calculus. The unnormalized average leads to a kind of fractional exclusion distributions. 
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